Finite element discretization and Newton-multigrid solution techniques for the log-conformation reformulation (LCR) of viscoelastic flow problems ⋆
نویسندگان
چکیده
We present special numerical techniques for viscoelastic fluid flow utilizing a fully coupled monolithic multigrid finite element approach with consistent edge-oriented stabilization technique. The governing equations arise from the Navier-Stokes equations for the Oldroyd-B type of fluid with the help of the log-conformation reformulation (”Kupferman’s trick”) to allow a wide range of Weissenberg numbers. The resulting nonlinear system consists of 6 variables for velocity, pressure and the logarithm of the conformation stress tensor in 2D. The system is discretized in time by using a fully implicit second order accurate time integrator. In each time step, we have to solve a discretized system in space employing the high order finite element pair Q2/P1. We utilize the discrete damped Newton method with divided differences for handling the Jacobian, and apply a geometrical multigrid solver with a special Vanka smoother to handle the linear subproblems. Local refinement can be assigned at regions of interest to reduce the computational efforts. The presented methodology is implemented and validated for several well-known benchmark problems.
منابع مشابه
FEM techniques for the LCR reformulation of viscoelastic flow problems
We present special numerical techniques for viscoelastic fluid flow utilizing a fully coupled monolithic multigrid finite element approach with consistent edge-oriented stabilization technique. The governing equations arise from the Navier-Stokes for the Oldroyd-B type of fluid with the help of the log-conformation reformulation to allow a wide range of Weissenberg numbers. The resulting nonlin...
متن کاملEfficient FEM Solvers for Incompressible Nonlinear Flow Models
We present special numerical techniques for non-Newtonian flows including non-isothermal, shear/pressure dependent and particularly viscoelastic effects which are based on the logconformation reformulation (LCR), including Oldroyd-B and Giesekus-type fluids as prototypical models. We utilize a fully coupled monolithic finite element approach which treats the velocity, pressure, temperature and ...
متن کاملMonolithic Newton-multigrid solution techniques for incompressible nonlinear flow models ⋆
We present special Newton-multigrid techniques for stationary incompressible nonlinear flow models discretised by the high order LBB-stable Q2P1 element pair. We treat the resulting nonlinear and the corresponding linear discrete systems by a fully coupled monolithic approach to maintain high accuracy and robustness, particularly with respect to different rheological behaviour but also regardin...
متن کاملNumerical simulation of a rising bubble in viscoelastic fluids
In this paper we discuss simulation techniques for a rising bubble in viscoelastic fluids via numerical methods based on high order FEM. A level set approach based on the work in [9] is used for interface tracking between the bubble and the surrounding fluid. The two matters obey the Newtonian and the Oldroyd-B constitutive law in the case of a viscoelastic fluid while the flow model is given b...
متن کاملA Newton Based Fluid–Structure Interaction Solver with Algebraic Multigrid Methods on Hybrid Meshes
Fluid–structure interaction problems arise in many application fields such as flows around elastic structures or blood flow problems in arteries. One method for solving such a problem is based on a reduction to an equation on the interface, involving the so-called Steklov–Poincaré operators. This interface equation is solved by a Newton iteration for which directional derivatives with respect t...
متن کامل